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MSTRACT 

In this paper we explore the consequences of particular stage linkage 
struciurea for the evolution of a population. We first argue the importance 
of constructing ''dynamic" models of developmental theories and show through 
a series of eKamples the implications of various stage connections for popu- 
lation movements* In discussing dynamic models, one thrust of our comments 
is to identify the sorts of process features about which assumptions must 
be made in order to convert a static theory about stage connections (the 
sort of specification conmonly presented in life-span psychology) into a 
dyiiamic model, A second focus of our discussiDn concerns inverse problems^ 
how to utilize a model formulation so that the stage linkage structure may be 
recovered from survey data of the kind collected by developmental psychologists. 
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MATHEMATICAL REPRESENTATIONS OF DEVELOPMENT THEORIES 

Burton Singer Seymour Spilerman 

Columbia University University of Wisconsin 

I, INTRODUCTION - 

Although tlmaj usually in the guise of age^ is a crucial variable in 
developmental psychology, it is the case that formal models of developmental 
phenomena rarely have the character of dynamic representations , in the sense 
of mimicking the evolution of an empirical process through tiiaey Th& analytic 
procedures employed most extensively by life-span psychologists are factor 
analysis s regression, analysis of variance, scallngs clustering, and variants 
of these methods (see, for Instance , Nesselroade and Reese 1973). These are 
powerful cechrtiques for identifying variables that are central to the course 
of development in a particular substantive area (e*g., intellectual matura^- 
tioni acquiiltion of moral values). AlsOj when applied to panel data, the 
procedures can yield insights into how the salience of key variables shifts 
over the life cyclej or over a portion thereof (e,g,j stages in infancyj youths 
adulthood). 

These analytic methods do not j however, lead to dynamic formulations of 
developmental theories^ which can be useful in tasting predictions from a 
theory about the evolution of an empirical process, or in comparing the 
implications of competing explanations* By a dynamic formulation we mean a 
representation which incorporates into the mathematics the main assumptions 
about a deveiopmental phenomenon and is specified in such a way that the 
relevant variables, and their postulated interrelations, are functions of 
.. time or subject's age. In this sense, like the empirical process, it too 
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constitutes an evolving system. As a simple Illustration of such a rnQdel, 
considar the following statements of alternative evolutionary mechanisms: 

(A) The growth of a process at each Instant is proportional to its 
potential for future growth. 

(B) Itie growth of a process at each instant Is proportional to Che 
product of itB nurrent size and its potential for future growth. 

These statements might be proposed as competing explanations of the 
Qianner by which information is diffused in a population of size N. In formu- 
lation (A), it matters not how many persons y(t) know the information of 
concern at instant t; only those yet to hear, numbering N-y(t), are salient 
to the diffusion rate. If the information were propagated by a mass media 
/source, such as radio or television ^ rather than by interpersonal communica- 
tion, this model might apply. Formulation (B) ^ in comparison, is consistent 
with a process' in which those already aware of the inforniation "infect" the 
uninitiated through contact and conversation. Assuming that the informed 
and the uninformed mix randomlyj the variable governing the evolution of the 



process would be y(t)[^N--y (t)J , which measures the rate at which individuals 
from the two groups come into contact. 

The evolutionary mechanisms^ (A) and (B) , can be represented by the 
differencial equati ons (1) and (2), respectively, 



where and are constants which adjust for the time unit (e.g, , day, year) 
used in the moasurementa . ^ Equations (1) and (2) have for solutions (3) and (4) 




(1) 



d(t) 




y(0) = 1 



(2) 



y(t) 



N<1 - a 




) 



(3) 



Ne 




y(t) 




(4) 



which predict the different evolutionary paths displayed in Figura 1. 

Thase formulations are "dynamic" in that time appears eKplicitly as a 
variable; they are process "models" in that the predicted value of y(t) evolves 
according to rhe assumptions of a particular theory. If a researcher has data 



(4)s or a specification of an equivalent sort best approKimates his observa- 
tions* By this exercise it is often possible to select among competing 
explanations of the mechanism underlying a developmental process* Indeed, 



data on drug adoptions by physicians (also see Coleman 1964, pp* 43-45). They 
concluded that the drug acquisitipn pattern by socially integrated MD*s is 
best represented by a logistic curve (implying mechanism [B]), while isolated 
MD's adopt according to the constant source model (mechanism [A]), as they are 
Influenced principally by drug advertisements In trade journals* To our 
knowledge, although developmental psychdlogists emphasize ontogenetic processes 
and employ the imagery of an evolutionary system, few attempts have been made 
to translate their theories into formal models of the above sort* 

In this paper, we describe the formulation of dynamic models where the 
objective Is to test developmental theories against data or ascertain the con- 
sequences of particular aasumptlons about the etructure of a process* To 



on the time course of an empirical process, he could test whether equation (3), 



these very models have been applied by Coleman, Katz, and Menzel (1957) to 
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Figure 1* Illustrative Growth Curves for Diffusion via Social Interaction 
and Diffusion from a Constant Source- 



y{t) 




N = population size; y(t) = number aware of the information 
at time t. 



ERIC 



8 



delimit our task, we focus on the sort of macheniatlcs that is appropriate 
for studying qualitative change. As a result, the tools we introduce are 
pertinent to thaorias which pogtulata stage sequences, a variety of aKplana= 
tion with considerable precedent in developmental psychology (Piaget 1960; 
Kohlberg 1968; Ausubel and Sullivan 1970)* To the degree possible we have 
organized this paper with a view toward substentive issues and have concen^ 
trated on the translation of theoretical specifications into mathematical 
formalism; the reader usually is referred elsewhere for mathematical details 
and estimation procedures. The organization of the paper is as follows ^ In 
the next section we introduce a class of models that is suitable for studying 
evolutionary processes that incorporate the notion of stage. In section III 
we describe how particular stage theories can be cast in the framework of the 
general model. In section IV we relax several requirements of the basic model 
so that it can more realistically represent developmental phenomena, 

II, THE CONCEPT OF DEVELOPtffiNT STAGES AND A MATHMATICAL FORMULATION OF 
STAGE PROGRESSIONS 

Stage sequences have been postulated for a variety of developmental 
processeS"~the evolution of moral behavior (Kohlberg 1973) , cognition (Piaget 
1954), peraQnality (Loevinger 1966), and motor skills (Shirley 1933), to cite 
but a few topics. There also exist diverse formulations of stage models in 
the literature of life-span psychology. These differ with respect to the 
presumed sources of the atagas and with regard to the rules governing movement 
between them. In regard to stage origins i some authors have emphasized 
maturatlonal considerations, in which individuals are viewed as programmed 
genetically for particular behaviors or abilities to emerge (Gesell 1954). 
The specification of paychosexual stages, keyed to biological activation of 
the seK glands j provides an Illustration (Kohlberg 1973, p. 181), Others 



view etagea as arising from intaractions with the social environment. Kohlberg 
(1968, pp. 1016-1024), for example, contends that experience with the cultural 
and physical world is necessary for cognitive atagas to taka the shapes they 
do. Still other rBsearchers have adopted the position that stages are a useful 
research construct around which to discuss development, without insisting that 
they have an empirical existence (Kaplan 1966| Reese 1970). 

We shall not discuss further the very important Issues concerning the 
etiology of stages, but will focus instead on the mathematical r«prasentatlon 
of theories about stage connections and on the consequences of various linkage 
structures for the evolution of individuals among the stages. Formulations 
of stage connections In a developmental process differ according to whether 
the progression Is viewed as unilineal or multilinaal, whether stages in the 
sequence can be skipped, and whether regression to an earlier level is possible. 
A second set of considerations pertinent to the structure of developmental 
theories concerns the age specificity of a stage and the related matter of the 
virlabllity of duration In a stage. For discussions of these topics in the 
context of particular substantive processes, the reader is referred to Emmerich 
(1968) and Kessen (1962). 

To develop the mathematical apparatus for ascertaining the Implications 
of particular stage connections, we discuss both the simplest prototype of a 
stage theory (for concretenass) and the general mathematical formulation. ^ ' 
Consider, then, a developmental progression consisting of n stages, in which 
r.he linkage Is unilineal and there is no possibility of stage skipping or 
regression. An example of such a structure, with n equal to 5, Is presented 
in panel A of Figure 2; henceforth this model is referred to aa example 1 . 
It will be convenient to also have available a matrix representation of the 
stage linkages. For an arbitrary n-stage structure, m define a matrix M, 

er|c 10 
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"^11 ^12 
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whose elements are m - {probability of transferring fram staie 1 to stage J 

when a transition OGCurs}, where 0 £m <_ 1, and E m ^ 1, These restrict 

IJ j^l 

tlons on the elements of M ensure that each row of the matrix constitutes a 
probability distribution. We require, in addition, that m^^ - 0 for each stage 
i which is not an absorbing state of the process; that is, from which individual 
can exit. This means we exclude the possibility of withln'Stage transitions , 
a type of move which is undefined in most developmental theories. Also, we set 
^ii ^ ^ stage which is an absorbing state of the process. This is 

done for mathematical convenience and, as we shall sea, carries no substantive 
implications. In the particular case of the unilineal progression (Figure 2, 
panel A), we have the further requirements on ^^^^^2. ^ l^and m^^ ^ 0 other- 
wisa (except that - 1). This matriK, M^, is reported in panel B of Figure 2 

To this point, though matrix M conveys important s actural information 
about the process, the description of the stage progression is a static repre- 
sentation. To elaborate the model we must indicate how stage transition events 
occur. At a general level of description we assume that the time spent by 
an Individual In stage i follows some probability distribution , 



Prob^(tj^ < t|tj^,..,.tj^^^) 



(6) 



where t ,...,t report the sojourn times in earlier stages. Our imagery, 
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Figure 2. Rfipresentatlon of a Simple Unlllneal Stage Structure 



a. Diagram of Stage Linkages 



© -^®— — 1^0- 



Matrix Representation* of the Stage Llnkag, 



es 



0 1 0 0 0 

0 0 10 0 

0 0 0 1 0 

0 0 0 0 1 

0 0 0 0 1 



Each row of K^^ Is a vector of destination probabilities. Thus, 
if an individual were In stage one before a transition, the row one 
entries would pertain and they indicate movement to stage two with 
probability equal to 1. " ' 

b_. . 

The main diagonal entries are set equal to zero (with the 
exception of row 5) to Indicate that a "move" Is not defined apart 
from a stage transition; i.e., there is no notion of movement within 
a stage. The main diagonal entry of row 5 is set equal to 1 because 
this stage la an absorbing state (m^^^ - 0 for J # 5) and the definition 
of M^— see text — requires £m - 1. 

j ^ 
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therefore, Is the following. An individual originates in stage 1 at 

the beginning of the process, t^ « 0, He remains there for an interval 

T^^ specified by a distribution function Prob^(T^ < t), and then transfers to 

stage j with probability * He remains in this stage for a period t^i speci^ 

fied by a conditional probability distribution Prob^ (tg < t|T^), then transfers 

3 

to stage k with probability m * and so forthi The process continues until somei 
absorbing state is reached, at which point the evolution la terminated. The time 
path for the unilineal progression associated with the stage linkages of matrix 
is presented in Figure 3* 

Several further assumptions are necessary to complete the specification 
of the model. One matter concerns the relevance of an Individual's past move-- 
mant history to the course of his subsequent evolution among the stages* We 
assinne, 

(1) knci^Sdge of current stage conveys all information that is 

relevant to forecasting future movements. 

Stated technically i if m . ^ - {probability of moving from stage 1 to 

ij J ab t • * t 

stage j at the occurrence of a transition, given prior sojourns in stages a, 
b, . , * f } , then 

"lj.ab...f " °1J • ■ ^ 

(This assumption is superfluous in the current example of a unilineal pro- 
gression since there Is only one possible path, but It Is relevant to the 
evolution of a population in less restrictive models.) We indicate in the 
next section that this specification has been employed in descriptions of stage 
linkages in developmental psychology* 
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Figure 3. A Sample Path Descrlptlou Corresponding to the Unllineal Staee 
Structure of Figure 2. 



Stage 



3 • 



2 ■ 



1-r 

I 



■^ 't^ 



time 



It is assumed that there are five stages, which must be traversed 

sequentially, la the value of a random variable and denotes the 

sojourn time for an Individual in stage i. Stage 5 Is an absorbing 
state of the process. 
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For an initial baseline class of models, we further asBuine, 

(ii) the sojourn time in stage i is Independent of previous sojourn 
times and is exponentially distributed; that is, 

-X t 

^"^l^^k * ^1^1"" 'Vl^ " Prob^(Tj_ < t) s F^Ct) = 1-e ^ . (7) 

Use of the exponential distribution amounts to specifying that the probability 
of departing from stage i during the infinitesimal interval t + dt, condi- 
tional on being In stage i at time t, equals 



-X t 

f^(t)dt X^a ^ dt 



X.dt 



1 - FiCt) -X.t " 

l-(l=e i ) 

where f^(t) ia tha density function corresponding to F^(t), This result, in 
tum^ indicates that the probability of leaving stage i is independrnt of 
duration in the stage, and is tantamount to specifying an absence of aging, so 
new entrants have the same livelihood of departing as individuals who have been 
in the stage for some period of time. The parameter, X^, incidentally, has 
an interpretation as the rate of movement out of stage i; consequently, 1/X^ 
equals the expected duration in stage i. 
Finally, we require that 

(iii) If the data pertain to the movements of a population, rather 
than to the transitions of a single individual, the popula- 
tion Is homogeneous with respect to the structure of the 
evolutionary process. 
This does not mean that all persons have the same duration in stage i, but 

that t^^, the time spent In stage 1 by individual c, follows the single expo- 

-X c 

nentlal distribution F^(t) » 1-e - . Stated less formnlly, duration in a stage 
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Is a random variable with the underlying distribution of holding times the same 
for all indlvidualfl. Similarly, where alternative destinations are available 
to persons in stage 1, homogeneity means that all have the same list of prob- 
abilities for making the various transitions, not that they move Identically. 

It Is worth dwelling on the conceptual statuB of the preceding assump- 
tions. The question of the structure of M Is a familiar topic to developmental 
psychologists, since stage theories are commonly specified at this level. 
Assumptions (1) to (iii) can be viewed as "side conditions," aspects of the 
process to which researchers have generally not been sensitive, though see 
Kessen (1970) and Eminerlch (1968) for provocative comments on precisely these 
matters. What Is made evident by formulating a dynamic model Is that develop- 
ment theorists must address these auxiliary questions if complete models are 
to be specified. The particular assumptions we have made constitute a gross 
simplification of reality; this is especially true of specification Xli) . which 
postulates an absence of duration effects, and specification (ill), which 
postulates population homogeneity. Wlese assumptions do, however, provide a 
convenient starting point from which to consider more realistic formulations, 
which are developed in the next sections. 

We now wish to convey the implications of assumptions (i) to (111) for 
the movements of Individuals among the stages. We denote by P^j (t) the prob- 
ability that an individual in stage i at time 0 moves to stage j by time t. 
(This probability differs from m^^ in that the latter refers to movement pro- 
clivities at the occurrence of a transition, not over widely spaced time 
intervals.) With this specification In hand, the evolution of a population 
among the stages Is described by the ayatom of integral equations, 

16 0 1 i. J 1 n 
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where 5^^ ^ 1 If i ^ j , and 0 otherwise. This eKpresslon, known as the back- 
ward equations for a CDntinuous-tima Markov process (Feller 1971, p, 484), is 
amenable to the following Interpretations (1) When i # j, p^j(t) consists of 
the sum of products of three factors! the probability of a ^Irst departure 
from stage 1 at time the probability of a stage i to stage k transition at 
that Instant, and the probability of transferring to stage j by some combination 
of moves in the interval t - u. The summation is over all intermediate stages 
k and over all time divisions u in the interval (0, t)* (2) When 1 ^ j, in 
addition to the above term, there is the possibility of not transferring out 
of stage i during (0, t)* This probability is given by the first term* 
If we represent by P(t) the matrix of elements p^jCt), 



nn 



0 £ p (t) £ 1, Zp,,Ct) ~ 1, then the integral equations (8) have the con- 
venlent solutionj 



P(t) 



A[M - Ijt 



P(0) - I. 



(9) 



In this representation A Is a diagonal matrix, 



ERIC 



^ J 

n 



whose encrles are the reciprocala of the expected duration times in each ataga^ 
I la the identity matrix, and M is the array Bpeclfled in equation (5) which 
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describes the pattern o£ movement between the scages. Further, by the expression 
A 

e , A an arbitrary squa^re matrix, we mean the power series In A, 



I 



n«0 



nl 



(10) 



which can be evaluated by standard numerical methods (see, e.g., Gantmacher 
tl960]). 

it is useful to recapitulate what is accomplished by this mathematical 
formulation. The matrix P(t) relates the distribution of a population among stages 
at time t to its distribution at time 0, In the sense that a typical entry. 
P^j(t), represents the probability of moving from stage 1 to stage j during the 
Interval (O.t). The model is "dynamic" in that P(t) is a function of time; 
with the passage of time P(t) describes the evolution of the population among 
the stages. Equation (9) shows how the matrix P(t) is built up from the arrays 
M and A. However, while this equation is useful as a calculating formula, the 
logic of the process is conveyed more adequately by the integral equations (8). 

To Illustrate this model in the setting of a simple unilineal progression 
(matrix M^^ of Figure 2). we must specify average waiting times in stages 1.2,3. 
and We assume these to be .5, 1, 2, and 5 years, respectively. Consequently, 
we have for matrix A, 



2 
0 
0 
0 
0 



0 

1 

0 
0 
0 



0 

0 
.5 
0 
0 



0 
0 
0 
.2 
0 



0 
0 
0 
0 



(11) 



where the choice of is arbitrary. (Since stage 5 is an abaorblng state, 
the notion of waiting time to a departure hns no meaning. Mathematically, 

- IJjj - tmjj - 1] - [1-1] - 0, BO bears no Influence on the calculations.) 
Now, from M, , A, and I, we have 
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A<Mj^ - I) 



-2 


2 


0 


0 


0 


0 


-1 


1 


0 


0 


0 


0 


-.5 


.5 


0 


0 


0 


0 


-.2 


.2 


0 


0 


0 


0 


0 



(12) 



For the illustrative times t « 1, 2, and 4 years, we obtain, from (9), for P(t)j 



PCD 



,1353 
.0000 
,0000 
,0000 
,0000 



.4651 
.3679 
.0000 
.0000 
.0000 



,3263 
,4773 
,6065 
,0000 
,0000 



.0691 
.1438 
.3537 
.8187 
.0000 



.0041 
.0110 
.0398 
.1813 
1.0000 J 



(13) 



P(2) 



,0183 
.0000 
.0000 
.0000 
,0000 



.2340 
.1353 
.0000 
.0000 

.0000 



.4641 
,4651 
,3679 
,0000 
,0000 



,2482 
,3394 
,5041 
,6703 
,0000 



.0354 
.0602 
.1281 
.3297 
1.0000 J 



(14) 



and 



P(4) 



.0003 
,0000 
.0000 
.0000 
,0000 



.0360 
,0183 
,0000 
,0000 
,0000 



,2881 
.2340 
.1353 
,0000 
,0000 



.4843 
.5079 
.5233 
.4493 
.0000 



.1913 
.2398 
.3413 
.5507 
1.0000 



These values of PCt) describe the evolution of individuals among the stages, 

subject to the assumptions about the process structure detailed above. The 

4 

entries of p^j(t) refer to proportions of the population who have moved 
between particular stages in the relevant time interval. For example, according 
to the entries In the top row of PCI), if observations are taken one year apart, 
we would expect l"* percent of the population in stage 1 at time 0 to still be 
there, 46 percent to have moved to stage 2, and 33 percent__to have reached 
stage 3. By comparison, over a four-year interval, less than 1 percent would 
reidfiln in Btage 1, 48 percent would have reached stage 4, and 19 percent would 
be in the terminal stage of the process, 

19 



16, 

The results from the three calculations reveal that, even though the 
progression is unillneal with all individuals characterized by the same para- 
meters, if observations were taken on the population at two time points, t^O 
and t^t^, the array p(t^) might be interpreted as evidence for a more complex 
theory, such as one permitting stage skipping or population heterogeneity in the 
rate or pattern of movement. Further, the correipondence between the matrix 
constructed from the population locations at two time points, PCtj^), and the 
rule governing stage transitions, M^^, decreases with time. Thus, different 
researchers observing the same population at two time points, but with different 
spacing intervals, might draw contrary conclusions about the stage linkage 
structure even though -the single mechanism, of Figure 2, governs its 
evolution. Only with a formal model of the process could one hope to uncover 
its underlying structure* 

III. MODELS OF MORE ELABORATE STAGE THEORIES 

The matrix M contains structural Information about stage linkages. Since 
theories of development are commonly posed at the level of specifying this 
array, fleKibillty in Incorporating a variety of specific formulatioris would 
appear to be an important feature of a general framework for describing 
evolutionary behavlorp In this section we focus on the issue of translating 
stage theories into M-matrlces, and illustrate the evolution of PCt), the 
transition matrlK for a population based on Its locations at times 0 and t, 
under altarnative specifications of As we have noted, auxiliary information 
about the process i concerning the distribution of waiting time intervals and 
the form of population hetProEenelty , is required for a full description of 
a dynamic modeli In the next soctlon we therefore elnborate upon these "side 
conditions" and outline ways In which our initial asaumptions can be relaKedt 
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No technical difficulties arise in reforaulatlng the continuous-time Markov 



model to accommodate more elaborate theories of stage linkages than the struc- 
ture In Figure 2, We Illustrate the procedure with a few examples,^ 

(2) A unlllneal progression which permits stage skipping . The formula- 
tion of such a structure is diagrammed in Figure 4, panel A; its translation 
into an M-matrlx Is reported in panel B. The principal new feature is that, 
supplementing the deterministic sequence of Figure 2, it is now possible to 
move directly from stage 2 to stage 4 and from stage 3 to stage 5, when a 
transition out of the relevant origin locntlon takes place. We must also specify 
the probabilities of following the alternate paths. In the present example, 
lacking information as to the relative magnitudes of the various probabilities , 
we assume all destinations to be equally likely; that is, we prescribe m^^ ~ 
m^^ " .5, and m^^ = m^^ « .5. In practice, estimates of the transition prob- 
abilities would be assigned on the basis of theory or from observation on the 
en^lrlcal process. 

Using matrix M^, together with the A array of equation (11), whose entries 
describe the rate of movement by individuals out of each stage, we obtain for 
and P(4), from equation (9), 



.1353 
.0000 
.0000 
.0000 
.0000 



.4651 
.3679 
.0000 
.0000 
.0000 



.1632 
.2387 
.6065 
.0000 
.0000 



.2012 
.3177 
.1768 
.8187 
.0000 



.0352 
.0757 
.2166 
.1813 
1.0000 



(16) 



P(4) 



.0003 
.0000 
.0000 
.0000 
.0000 



.0360 
.0183 
.0000 
.0000 
.0000 



.1440 
.1170 
.1353 
.0000 
.0000 



.4104 
.3964 
.2617 
.4493 
.0000 



.4093 
.4683 
.6030 
.5507 
1.0000 



(17) 



18 



Figure k. Representation of a Unilineal Progression in which Stage Skipping 
Is Permitted rr e 



a. Diagram of Stage Linkages 




b. Matrix Representation of the Stage Linkages* 



M. 



0 10 0 
0 0 .5 .5 
0 0 0 .5 



0 

0 



0 
0 



0 
0 



0 

0 



0 

0 
.5 
1 
1 



All destination stages corrosponding to an origin location are 
.8su.ed to occur with aqual probability. See notes to Figure 2 for 
additional details on Interpretation of Mg. 
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These P(t) arrays are the transition matrices a researcher should expect to 
observe if the stage locations of individuals are surveyed one year or four 
years apart, assuming that the population evolves according to the linkage 
specification together with the auxiliary conditions outlined in the pre- 
ceding section. The entries are different from those obtained with the simple 
unllineal progression (equations 13 and 15) ^ yet the same pattern of zero's 
and non-zero's is present ^ and without a formal model of the evolution of the 
process a researcher would be unable to predict the different implications of 
these structures. 

(3) A unllineal progression with stage skipping and the possibility of 

regression * We now superimpose on the linkage structure the possibility of 

reverting to an earlier stage. This arrangement is diagramned in Figure 5, 

panel in which we have /rovided for the possibility of bac^ard flows from 

stage 2 to stage 1, from stage 3 to stage 2, and from stage 5 to stage 4, The 

M-matriK corresponding to this model is reported in panel B* Again j where 

multiple destinations correspond to an origin stagey we have arbitrarily assigned 

eaual values to the m, .'a. There is one additional alteration in M > in com- 

Ij 3 

parison with the M-mat rices of earlier examples. Because there now exists 

a possibility of regressing from the terminal stage to an earlier level, m^j # 1. 

To maintain our conceptual imagery, In vhich within-stage transitions are 

undefined, we set m^, ^ 1 and m^^ ^ 0. Note that the former value does not 
' DH 55 

imply a high rate of departure from stagti 5, since the rate of movement is 
controlled by X^, It only means that all iransltions from stage 5 are directed 
to stage 4* 
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Figure 5. Representation of a Unillneal Progression In which Stage Sklppln 
and Regression to an Eairller Level are Permitted 



a. Diagram of Stage Linkages 




b. Matrix Representation of the Stage Linkag 

0 1 0 0 0 ' 

.33 0 .34 .33 0 

^2 ^ 0 .33 0 .34 .33 

0 0 0 0 1 

0 0 0 1 0 



All destination stages corresponding to an origin location 
are aasumed to occur with equal probability. See notes to Figure 2 for 
additional details on Interpretation of M . 
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To obtain P(t) we use and A in conjunction with equation (9). Here 
the element in equation (11) is no longer arbltfary, as movement ouw of 
stage 5~£s~"a possibility. We shall assune that such reversions are rare, and 
hence specify the average waiting time to a transition from stage 5 to be 
eight years; that is, « .125, With these assumptions, we obtain for our 
illustrative calculations at t = 1, 4, 



P(l) 



,2043 
,0871 
,0094 
,0000 
,0000 



.5240 
.4758 
.0858 
.0000 
.0000 



,1153 
,1742 
,6215 
,0000 
,0000 



.1374 
,2217 
,1461 
,8292 
,1067 



,0190 
,0411 
,1371 
,1708 
.8933 



(18) 



P<4) 



.0330 
.0259 
.0135 
.0000 
.0000 



,1560 
,1246 
,0739 
,0000 
,0000 



,1652 
1500 
,1846 
,0000 
,0000 



.4025 
.4174 
.3542 
.5523 
.2798 



,2433 
.2820 
,3738 
,4477 
.7202 J 



(19) 



If we compare the P(l) mattices and the P(4) matrices from the three 
examples [i.e., equations (13;, (16), (18) and (15), (17), and (19)], we can 
acquire a fair idea of the implications of different stage Interconnections 
for the evolution of a population among the statuses. We also emphasize the 
fact that if a population were surveyed at two time points, especially widely 
spaced time points, it may not be obvious from inspecting the empirically 
deterTDlned transition array, P(tj^), as to the structure of the stage linkages 
(matrix M) which generated the observacions. We will return to the issue of 
Identifying the correct structure and recovering matrix M when the observations 
on a process are widely spaced; first we conclude this discussion on translat- 
ing theoretical specifications of stage linkages into M-matrices with a couple 
of exnmplea of multillnonl sequences that have been described In the develop- 
mental psychology literature. 
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<*) A diverge nt multiple progression (Van Den Daele 1969, Figures 2, 4), 
nils stage linkage structure has the diagrammatic representation of Figure 6, 
panel A; its corresponding M^matrix is presented in panel B. Because stages 
4'-7 are specified to be terminal scates of the process, the corresponding rows 
of have I's in the main diagonal. Van den Daele provides no discussion 
of waiting time distributions to departure from the various stages; hence the 
model remains incomplete as an evolutionary process. 

(5) A convergent multiple progression (Van Den Daele 1969, Figure 2). 
This stage sequence is depicted In Figure 7, panel A, and its associated M- 
matrix is reported in panel B* In this Instance, the structure consists of a 
collection of deterministic unilineal progressions, the specific sequence for 
an individual being contingent upon his entry stage* Note also that the 
assumption of irrelevance of past hlstoryj which is posited in this formula^- 
tlon. Is one of the side conditions we have required (assumption [i] in the 
preceding section). In particular^ this specification appears in the fact that 
knowledge of the path by which one has reached stage 5 {^.^ stage 6) is of no 
value in forecasting^ or understanding, an individual's subsequent movementa* 
Van Den Daele (1969) discusses several additional models of stage linkages, 
such as ^'partially convergent, divergent progression,*' and "partially divergent * 
convergent prDgression. " As the procedure in converting flow structures Into 
M--matrlces should be evident at this point, discussions of these specifications 
are not presented. 

To recapitulate, subject to several side conditions, we have shown that 
it is possible to construct formulations of a range of developmental phenomena 
which mimic the evolutionary character of the observed process. With such a 
model one can forecast the movements of a population among the stages. By 
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Figure 6. Representation of a Divergent Multiple Progression' 



a. Diagram of Stage Linkages 




Matrix Representation of the Stage Linkages 
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^Source I Van Den Daele (1969, Figures 2, 3), 
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Figure 7. Represantation of a Convergent Multiple ProgreBslon^ 



a. Diagram of Stage Linkages 




h. Matrix Representation of the Stage Linkages 
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Source I Van Den Daele (1969 , Figure 2), 
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carrying out the requisite calculationa for different specifications of the 
stage linkages, and comparing the predictions, It Is possible to ascertain the 
ways in which rather complex theories produce divergent implications and design 
testing schemes which maximize the possibility of rejecting one or another 
formulation as a description of the empirical process. Of equal importance, 
it is often possible to work backwards , starting with observations on the stage 
locations of a population at a few widely spaced time points, and derive the 
structure of the stage linkages compatible with the data. 

An inverse problem . Until this point we have assumed that observations 
have been made on an empirical process in a way such that M and A can be esti- 
mated directly from the data, or that theories are available which specify 
the values of their entries, and have sought to derive the evolution of the 
process subject to the presumed structure. In developmental psychology, it 
is not uncomnon for a researcher to have many observations on a few individuals 
(e*g. ^ Fiaget 1954). Such a data collection scheme approximates "sample path 
information," a complete history on movements and waiting times of the sort 
Illustrated in Figure 3. Detailed observations on a few subjects is a research 
strategy not without its costs, however. One learns little about the frequency 
of rare events (e,g,, regression to an earlier stage, stage skipping, rare 
development paths) and acquires only the most rudimentary knowledge about the 
variation of duration times in a stage. It is therefore not surprising that 
investigators who rely on this approach tend to be oriented to uncovering 
universal rules (e,g,, Plaget 1960) rather than to elucidating individual 
differences and ascertaining the variety of developmental patterns. 
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Partly because of the llraitatlons of small data sets, it is becoming 

increasingly common to employ survey methods , In which a large population, 

sometimes thousands of individuals, is observed (or interrogated) at a very 

few time points (e.g., Baltes and Nesselroade 1972). The spacing intervals in 

.. such panel studies are usually , wide, often one or more years elapses between 

interviews, so it is not unusual for some subjects to have made multiple moves 

while others have made one or zero shifts between stages. The transit ion 

matrices which can be constructed directly from such observations are PCt)-- 

arrays, rather than M--arrays, and the stage linkages may not be readily dis- 

cernable. Indeed, determination of the movement structure which underlies the 

evolution of the population can be a difficult task. 

One approach to ascertaining the stage linkages from survey data involves 

consideration of the -'inverse problem" to the mathematical formulation of the 

evolutionary model (equation 9), Stated formally, we have available the 

matrix ?(t^) , constructed from observations on the stage locations of individuals 

at times 0 and t. , The typical entry in this raatrlK is p (t,) ^ n (t )/n whe 
-L ij 1 ij 1 i, ' 

n ^ {number of individuals in stage i at time 0} and n (t ) ^ {number of 
1* ij 1 

parsons who started in stage 1 at time 0 and are in stage j at time t^}* We 
wish to inquire whether it is possible to recover a unique M^matrix for the 
process and, where the answer Is af f irraative, we wish to estimate this matriK, 

inie first step in solving the Inverse problem is to take the logarithm 
of both sides of equation (9) , 

Q - A[M - I] - ^ ln?(t ) , (20) 

Just what we mean by the logarlthra of matrix P(t^), the conditions under which 
a solution to equation (20) will exist, and the circumstances under which the 
solution will be unique, are complex iaaues which are discussed at length in 
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Singer and Spilerman (1976) • Assuming we can obtain a valid and unique 
Q^matrlx from these calculations ^ a second task, separating M from A, still 
remains p In many instances, though^ this matter is of little concGrn, 
since the pattern of zeros and non^zeros in Q and M I will be Identical and 
deyelopraent theories are often posed_ at the level p| identifying permissible 
transitions. Moreover, because zeros are typically present in many main diagonal 
cells of M in models of developmental structures, a complete or near complete 
separation between M and A can frequently be effected. 

We conclude this section with an example of the calculations associated 
with the inverse problem. Suppose observations taken on a population at t ' es 
0 and have produced the transition matrix, 



.0224 


.2633 


.2402 


.1261 


.3479 


.0063 


.1758 


.2460 


.1735 


.3983 


.0216 


.0288 


.3758 


.5060 


.0679 


.0365 


.0745 


.0288 


.6794 


.1809 


.0005 


.0960 


.0460 


.0177 


.8397 



(21) 



Such data would appear to be consistent with a variety of evolutionary mechanisms* 
From inspection of PCt^^) we do know that regression to some earlier stage must 
be possible, otherwise all entries below the main diagonal would be zero. 
Little else about the structure of M, however, can ba inferred from inspection 
of P(tj^). Indeed, because of the sizable non-zero elements in most cells of 
the matrix, a researcher might conclude that direct transitions are possible 
between most pairs of stages. 

If we are willing to assume that matrix was generated by a continuous^ 

time Markov process; that is, via the evolution of the structure P(t) - e^^^ $ 
for some matrices A and M which satisfy the definitional restrictions enumerated 
in connection with equations (5) and (9), we can solve for A(M-I)tj^ using 



equation (20), This yields the array. 
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A(M - l)t, 





-4. 


4. 


0. 


0. 


0. 






0. 


-2. 


1. 


0. 


1. 






0. 


0. 


-1. 


1. 


0. 






0.20 


0, 


0. 


-0.40 


0.10 






0. 


0. 


25 0. 


0. 


-0.25 ^ 




aiid 


M can 


be 


separati 


id by 1 


smploying the 



(22) 



argument. From our earlier examples we know that a main diagonal element m^^ 
of M will equal zero If any off diagonal entry In the same row, m. , la differ- 
ent from zero. According to equation (22), each row of matrix M must have at 
least one non-zero off diagonal element; therefore m^^ = 0 for all values of 
1. With this Information we can obtain At^ uniquely. 



At. 



4 
0 
0 

0 
0 



0 
2 
0 
0 
0 



0 
0 
1 
0 
0 



0 
0 
0 

.4 
0 



0 
0 
0 
0 

,25 



(23) 



and soJviug for M provides the structure reported In Figure 8, panel A, 
The schematic representation of the stage linkages Implied by is shown in 
panel B, in which probabilities of the various moves have been appended to 
Che paths. 

The point to be emphasized Is that it is not apparent from Inspecting 
matrix P(t^) In equation (21) that the underlying stage linkages are those 
reported In Figure 8^ nor would any static analytic procedure be likely to 
lead a researcher to the correct conclusion. What is necessary Is to construct 
a model of the evolution of the process and solve the implied inverse problem 
for the parameters which correspond to the particular data set* (In the present 
aKample, we have assumed that the underlying model is a contlnuous-tlmo Markov 
process [i,e,, specifications (l)-(lll) of the preceding section] and have 



iolved for the matrices At^ and M which are compatible with the observed 
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^Tha process is assumed to evolve according to a continuous-time 
Markov f ormulatlDn- 

^Entries Indicate the probalklllty of a stage 1 to stage j move 
when a tranaltion takes place, 

^Probabilities of the various transitions are attached to the 
appropriate paths. 



30 

array P<t^). in that they would have given rlBe to this array if the postulated 
evolutionary proceaa were approximately correct, 

IV. ALTERNATIVE SPECIFICATION OP THE SIDE CONDITIONS 

In this section we discuss relaKing two of the more burdensome specifica- 
tions of the model, in the aanse chat they are likely to be inappropTtate as 
characterizations of developmental prooesses. We first consider the require- 
ment that the duration intervals in a stage must follow an exponential dis- 
tribution (assumption (ii) of section 11). Following these comments we turn 
to the requirement that the population be homogeneous with respect to the 
process parameters A and M (assumption (iii)). 

More general waiting times than ex ponential . The exponential distribution 
is frequently employed In the literature of reliability theory to describe 
duration intervals in a system state (stage In the current application). It 
has the advantages of being mathematically tractible and approximating reality 
in situations where the probability of a state change Is uninfluenced by 
aging or time in the state. For example, if the process states are "alive" 
and "not alive," then over the middle age ranges of many animal species, the 
age-specific mortality rate is relatively constant and the duration intervals 
(in the "alive" state) are reasonably well captured by the exponential dis- 
tribution. Similarly, when mortality results from exogeneous events— accidents— 
the distribution of ages at failure can often be approximated by the exponential. 

In a great many situations in social research, however, we know that 
pronenass to changing state is a function of duration. In particular, this 
has been suggested with respect to residence location (McGinnls, 1968) and 
eraployment affiliation (Ginsburg, 1971). In these applications It has been 
argued that the duration-specific departure rate decreases with time, giving 
rise to the phenomenon of "cumulative inertia"— the longer an individual 
remains in a state the less likely he is to leave in the inmediate future. 
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The eubstantive explanations for a declining daparture rate Involve the growing 
investment an individual has made, with duration, in friendships (in the first 
Instance) and in seniority in his place of work (in the second). There is 
no mathematical reason, however, to assume a declining departure rate In 
choosing F^Ct), and in other substantive contexts a different specification may 
be more appropriate* For a superb review of stochastic models incorporating 
the notion of duration dependence, see Hoem, 1972. 

A convenient way to generalize the ftorkov model to accommodate a variety 
of duration-time distributions is to begin with the Integral equation repre- 
sentation for transition probabilities* Equation (8) is a special case of 
the formulation 5 

t 

Pij(t) - fij_j[l - F^(t)] + I f^(u)m^j.pj^j(t-u)du 0 < l,j < n (23) 



in which the terms are identical with those of the earlier equation except that 

-X u 

f (u) replaces tha exponential densityp X e ^ and F (t) [the distribution 
function corresponding to f^Ct)], replaces [1 = e ]- A theoretically 



appropriate choice may now be made for F^(t), 



As an illustration, one candidate for F^Ct)^ in the case of a declining 
departure rate* is the two-^parameter family of functions 

-X t 

F^(t) « 1 - e - X^ > 0; 0 < < 1 . (24) 

Here the probability of departing from state i during the inf inlteslmfll interval 
(t, t + dt), conditional on the process being in state 1 at time t, equals 

Y.-i -x.t ^ 
f.<t) (X Y t )e . Y^-l 

1-F,(t) Yi 
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Because of the restriction on in equation (24), t * is a decreasing 
function of time, and the declining failure rate aspect of the distribution 
is evident. 

The general formulation (23) for duatlon time distributions and transi- 
tions between states generates a class of models known as semi-Markov processes 
These generally do not have simple repreaentations for the matrices P(t) 
analogous to aquation (9), and the solution of the systems of equations (23) 
requires numerical integration methods. 

Populatio n heteroRenelty . To this point we have assumed that the matrices 
A and M of equation (9) are identical for all individuals. This does not 
mean that all persons move identically since the process is probabilistic; it 
does imply, though, that Individual level characteristics are unrelated to the 
structural parameters of the process. In other words, homogeneity means that 
considerations of genetic makeup, intelligenca, sensory stimulation, and other 
factors by which Individuals differ from one another do not portend distinct 
evolutionary paths in the developmental process under consideration. 

There Is reason to believe, however, that individual differences are 
present in the course of davalopment in many processes (Werner, 1957; Kohlberg, 
1968, p. 1024), We therefore desire a formulation in which the movement 
pattern is parametrized in terms of variables which differentiate among persons. 
To construct a general specification of heterogeneity within the conceptual 
framework of a Markov process, we assume that, correspdndlng to equation (9), 
the atage transitions by individual c have the structure 

A (M - I)t 

^(t) - a ^ ^ . (25) 

This formula indicates that each person is characteriEed by a pair of matrices, 
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A and M , and his evolution, in turn, is describad by P (t). Thus* our 
formulation begins with a separate Markov process for each individual. 

This approach directs a researcher to identify the variables which describe 
heterogeneity; that is, to ascertain which factors account for individual 
differences in the matrices M and A, Thus^ not only does a heterogeneity 
formulation lead to more reaiistic models of evoiutionary processes, in that 
allowance is made for Individual differences , but it stresses the analytic 
tasks of specifying the variety of developraental patterns in a population^ 
and ascertaining the attributes which make an individual more prone to follow- 
ing one set of paths rather than another. 

One form of heterogeneity concerns the distribution of M=mat rices in a 
population. Focusing on these arrays serves to emphasise individual differ-' 
ences in proneness to making particular moves when a transition Cakes place • 
We shall not discuss this form of heterogeneity in the present essay and 
direct the interested reader instead to McFarland (1970) , Splleraan (1972a), 
and Singer and i Spilerman (1974)* A second form of heterogeneity stresses 
Individual differences in the A-matriXs i,e*3 in the rata:% <at Which departures 
occur for persons in the various states, We conclude this section with a simple 
formulation of population heterogeneity in which It is assumed that the indl= 
vidual differences can be expressed in the latter wayt 

To simplify the discussion i we further require the non--zero entries in 
the diagonal matrix A to be equal for an individual; ise*, ^ X for all 1# 
This means we are specifying Identical departure rates from all states. As a 
results equation (9) reduces to 0: 




e 



Xt(M-l) 



(26) 
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rate of movement value equal to A. We shall aasume that aquation (26) 
describes the evolution of an individual drawn at random from the popu- 



latlon. 



Heterosanalty Is Incorporated Into the formulation by apaclfying a density 
(-notion g(X) Which describes the distribution ot A-val„es In the population, 
le no„ define tha ioEHWlojrte^ -trl. corraspondlng to tl„as 0 



and t to be 



(27) 



^18 forn.ula expreesea the population-leval n^tri. as a weighted average of 
the lndividual=level arrays P(t|.). the weights reflecting the population 
Prcsu-Llona associated with particular X-values. 

to complete this specification of heterogeneity it is necessary to select 

.density function g(X) to describe the distribution of .-values. One useful 

choice is the gamma family of functions , 



„a a-1 -gx 
'(.-) ^ PA e 

r(a) ^ > 0, a > 0, B > 0 



(28) 



which is fl..ible enough to describe a variety of unimodal curves. With this 
selection of gCA), a convenient representation of the populatlon=level matrix 
P(t) is obtained (Spilerman, 1972b, p. 608)^ 

Il-e transition probabilities (29) do not describe the evolution of a fUrkov 
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process; however, they do describe the movements of a population In which each 
individual follnw^ ^ Markov model , with individual differences being specified 

by i(X) In equation (28). 

In analogy with our earlier Inverse problem discussion for Markov chains, 
the present formulation can be used with observations taken at widely spaced 
time points. 0 and t^;, togethex «ith e^imates^ of a and to yield an estlmat 
of the underlying transition mechanism M. according to the matrix equation 



„ „ ^-^-^l I fw. -^f''^ . (30) 



Thus, from observations of the sort collected In many surveys, even under an 
assumption of population heterogeneity in the rate of movement, it may be 
possible to recover the matrix of stage linkages which governs the evolution 
of the process. 

V. CONCLUSIONS 

In this paper we have explored the consequences of particular stage 
linkage structures for the evolution of a population. One thrust of our 
co^ents has been to identify the sorts of process features concerning which 
assumptions must be made in order to convert a static theory about stage 
connections into a dynamic model. A second focus In our discussion has 
centered on Inverse problems; how to utllizp. a model formulation so that the 
stage linkage structure (matrix M) may be recovered from survey data of the 
kind usually collected by developmental psychologists. 

We have presented only the most rudimentary sorts of stage structures, 
indeed, even within the Markov framework we have limited our consideration to 
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a subset of these models; namely, those which are time-stationary (l.e.i A 

and M are not functions of time). By this specification we have excluded 

the possibility of accommodating age-depandent transition lawSj a considera-- 

tion of substantial importance In devalopmental psychology, (An eKtension 

of the models discussed here to Incorporate both age dependence and cohort 

effects Is, however, a feasible imdertaking but with an increase In mat he--' ^^^^ 

matlcal complexity.) Further, all the models we have diacuesed entail a 

low dependence of future movements on the transition history of an individual ^ 

given his current stage. ^ Restrictions of these sorts are likely to be 

reasonable for some processes, unreasonable for others. Appropriate models 

of developinental phenomena must therefore be nf' ^f^-,ructed from a list of 

known characteristics about an empirical procesL 

We also point out that the concept of stage mnrges with the notion of 

state as the number and sorts of permissible transitions is increased, "Stage" 

seems conceptually rooted to the idea of progress (l*e, , development) and 

would be an appropriate component of a theory which sees the system's statuses 

as genetically determined or as facilitating the conditions for succeeding statu 

8 

to come into play. The mathematical framework we have introduced is also 
compatible with a "state" notion^ in which there is an extensive opportunity 
to cycle among the statuses* State formulations have been suggested in the 
psychology literature In relation to anxietys moods, etc. (e.g.^ Kessen, 
1962, pp. 72=^73). 

As a final set of considerations in relation to the structure of stage 
models, we note that all the formulations we have addressed are models of 
solitary processes. We have proceeded as if intelligence ^ cognition, motor 
skills, and personality development unfold autonomously. In reality there 
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no doubt exist eKtensive dependencies among some of these processes. Mathe- 
matical models of interacting developmental phenomena could be formulated 
but clear empirlcally-based specifications of such dependencies are still 
lacking^ 



4- , 
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NOTES 

^The Initial condition, y(0) - 1, In equation (2) is necessary because 
diffusion through communication cannot begin until at least one person Is 
knowledgeable. 

For a more technical presentation of continuous-time Markov processes 
see Feller (1968, Chap, 17) and Singer and Spilerman (1974). For discussions 
on the superimposition of theoretical structures on stochastic models see 
Coleman (1964, Chaps. 5, 6). 

3 

In the present example i, J,k=l. 2, 3, respectively. 

4 

If the observations are on a single individual the Interpretation of 
P^jCt) Is In terms of the probability of a stage 1 to stage j move between 
tlmas 0 and t. 

^The symbol over a matrix or over an element in a matrix, will mean 
that it should be viewed as estimated directly from data rather than calcu- 
lated from a mathematical model, 

^We begin here with example 2. Example 1 refers to the structure in 
Figure 2, 

^The time-stationary Markov formulations postulate an irrelevance of 
prior stage affiliations, durations In those stages, and duration in current 
stage. The last two of those restrictions can be eliminated by Introducing 
non-stationary semi-Markov models as deacrlbed, for example, in Hoem (1972). 
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Stagea In childhood, such as "walking'^ or "reading," expose an Indi- 
vidual to entirely new seta of eKperlences which may be prerequisites for 
the onset of more advanced behaviors. 
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